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Highlights 


We  present  several  new  formulations  for  the  classical  Capacitated  Vehicle 
Routing  Problem. 


We  survey  other  previous  flow-based  formulations. 

We  show  that  our  new  formulations  have  desired  properties. 


A 


We  prove  that  the  linear  programming  relaxations  of  the  n effi  formulations 
give  better  bounds. 
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Abstract 

The  Capacitated  Vehicle  Routing  Problem  is  a much-studied  (and  strongly  MV- 

\ 

hard)  combinatorial  optimization  problem,  for  which  many  integer  programming 
formulations  have  been  proposed.  We  present  two  new  multi- commodity  flow 
(MCF)  formulations,  and  show  that  they  domifipfe  all  of  the  existing  ones,  in  the 
sense  that  their  continuous  relaxations  yiel<pstronger  lower  bounds.  Moreover, 
we  show  that  the  relaxations  can  be  strengthened,  in  pseudo-polynomial  time,  in 
such  a way  that  all  of  the  so-called  knapsack  large  multistar  (KLM)  inequalities 
are  satisfied.  The  onlwother  relaxation  known  to  satisfy  the  KLM  inequalities, 
based  on  set  partitioning,  is  Strongly  A/T’-hard  to  solve.  Computational  results 
demonstrate  that  the  new  MCF  relaxations  are  significantly  stronger  than  the 


sea 


Vehicle  Routing  Problems  (VRPs)  are  classic  problems  in  operational  re- 
search and  logistics,  and  have  also  received  a great  deal  of  attention  from  the 
combinatorial  optimization  community.  A huge  number  of  papers  have  been 
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written  on  the  theory  and  applications  of  VRPs,  and  on  exact  and  heuristic 
solution  methods  for  them  (see,  e.g.,  the  edited  volumes  [5,  14,  27].) 

This  paper  is  concerned  with  the  Capacitated  VRP  (CVRP),  which  Dantzig 
and  Ramser  [7]  defined  as  follows.  A fleet  of  identical  vehicles,  with  limited 
capacity,  is  located  at  a depot.  There  are  n customers  that  require 
Each  customer  has  a known  demand.  The  cost  of  travel  between  any>< 
customers,  or  between  any  customer  and  the  depot,  is  also  known, 
find  a minimum-cost  collection  of  vehicle  routes,  each  starting  and  ending  at  the 
depot,  such  that  each  customer  is  visited  by  exactly  one  vehicle,  and  fio  vehicle 
visits  a set  of  customers  whose  total  demand  exceeds  the  vehicle  capacity. 

Letchford  & Salazar  [22]  surveyed  and  compared  several,sMteger  program- 
ming formulations  of  the  CVRP.  These  include,dLthe  so-called  two-  and  three- 
index  formulations,  the  single-,  two-  and  multi- commodity  flow  formulations, 
and  the  set  partitioning  formulations.  At  present,  the  most  successful  exact  al- 
gorithms for  the  CVRP  are  based  on  the  two-index  formulation  (e.g.,  Lysgaard 
et  al.  [23])  or  on  set  partitioning  formulations  (e.g.,  Fukasawa  et  al.  [10],  Baldacci 
et  al.  [3]). 

One  way  to  measure  tHe  s'ttength  of  an  alternative  formulation  is  to  project 
the  feasible  region  o^fts  continuous  relaxation  into  the  space  of  the  natural 
(two-index)  formulation.  Gouveia  [15]  showed  that,  in  the  case  of  the  single- 
commodity flow  formulation,  the  projection  satisfies  a family  of  valid  inequal- 
ities now  known  as  generalized  large  multistar  (GLM)  inequalities.  Letchford 
& Salazat  [22]  showed  that  the  projection  of  the  set  partitioning  formulation 
(with  only  elementary  routes  permitted)  satisfies  the  so-called  knapsack  large 
multistdr  (KLM)  inequalities,  defined  in  [21].  The  KLM  inequalities  include 
the  GLM  inequalities  and  the  so-called  subtour  elimination  (SE)  inequalities  as 
ecial  cases.  Unfortunately,  the  continuous  relaxation  of  the  set  partitioning 
formulation  is  itself  strongly  A/T’-hard  to  solve. 

This  paper  has  four  main  contributions.  First,  we  show  how  to  strengthen 
the  two  best  multi-commodity  flow  (MCF)  formulations,  by  adding  only  a poly- 
nomial number  of  additional  constraints.  Second,  we  show  that  the  projections 
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of  our  two  formulations  satisfy  the  GLM  and  SE  inequalities.  Third,  we  show 
that  the  new  formulations  can  be  further  strengthened,  in  pseudo-polynomial 
time,  in  such  a way  that  all  of  the  KLM  inequalities  are  satisfied.  (We  remark 
that  no  polynomial  or  pseudo-polynomial  time  separation  algorithm  is  known  fo: 
the  KLM  inequalities  themselves.)  Finally,  we  present  some  computational  re 
suits  that  demonstrate  that  the  new  MCF  formulations  are  significantly  stron; 
than  the  previously  known  ones. 

A,  mentioned  above,  the  current  best  algorithms  fo,  the  f*^based 
on  two-index  or  set  partitioning  formulations.  The  contribution  of  this  paper 
may  therefore  appear  to  be  only  of  theoretical  interest.  We  would  like  to  point 
out,  however,  that  there  exist  variants  of  the  CVRP  for  whi^JPit  is  natural,  or 
even  essential,  to  use  additional  commodity-flow  variables.  This  includes,  for 
example,  the  problem  described  in  in  [17],  if|  which  several  distinct  products 
have  to  be  picked  up  and  delivered  at  various  locations,  and  the  one  described 
in  [19],  in  which  the  cost  of  traversing;' an  arc  is  an  increasing  function  of  vehicle 
load.  Potentially,  our  results  could  be  used  to  derive  better  formulations  and 
algorithms  for  such  problems. 

The  structure  of  the  paper  is  as  follows.  The  literature  is  reviewed  in  Sec- 
tion 2.  The  strengthenea  MCF  formulations  are  presented  and  analysed  in 
Section  3.  The  result  on  KLM  inequalities  is  given  in  Section  4.  Some  compu- 
tational results  ate  given  in  Section  5,  and  some  concluding  remarks  are  made 
in  Section  6. 

Throughout  the  paper,  we  use  the  following  notation.  We  have  a complete 
directed  graph  G with  node  set  V = {0, 1, . . . , n}  and  arc  set  A.  Node  0 repre- 
sents thfe  depot,  and  nodes  1, . . . , n represent  customers.  We  sometimes  write 
Vc  for  V \ {0},  the  set  of  customer  nodes.  The  (positive  integer)  demand  of 
stomer  i e Vc  is  qt.  The  (positive  integer)  vehicle  capacity  is  Q.  The  (non- 
negative integer)  cost  of  traversing  arc  (i,  j)  G A is  Cij.  (Our  approach  can 
easily  be  adapted  to  the  case  of  symmetric  costs  and/or  the  case  in  which  the 
number  of  vehicles  is  restricted.) 
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2.  Literature  Review 

As  mentioned  above,  many  formulations  have  been  proposed  for  the  CVRP. 
For  brevity,  we  review  only  ones  of  relevance  here.  Subsections  2.1  to  2.4  cover 
two-index  vehicle  flow , single-  and  two-commodity  flow , multi- commodity  fL 
and  set  partitioning  formulations,  respectively. 

Laporte  & Nobert  [20]  presented  what  is  now  called  the  tj^dndex  vehicle 
flow  formulation.  For  all  ( i,j ) G A,  define  a binary  variabl^aq j , taking  the 
value  1 if  and  only  if  some  vehicle  travels  from  i to  j . For  any  S C V,  let 
S+(S)  (respectively,  6~(S))  denote  the  set  of  arcs  (i,j)  wit^f  G S,  j G V \ S 
(respectively,  with  i £ V \ S,j  €.  S).  If  S'  then  we  will  write  5+(i) 

and  S~(i)  rather  than  <5+({*})  and  4 ({/}}.  for  breviiy.  Given  some  F C A 


let  x(F)  denote  j)eFxij-  Finally, 
Q(S)  = Sies  qi.  Then  the  formulati^ft  is 


*))  = ! 
> \q(S)/Q] 

e {0,1} 


S6t  of  customers  S C Vc,  let 


(1) 

(2) 

(3) 

(4) 

(5) 


equations  (2)  and  the  in-deqree  equations  (3)  ensure  that  vertices 
a^Vi^dJ^actly  once.  The  constraints  (4),  called  rounded  capacity  (RC) 
equalities,  prevent  the  existence  of  infeasible  routes,  and  also  have  the  side- 
effect  of  preventing  subtours.  Finally,  (5)  are  the  integrality  conditions  on  the 
-variables. 

Several  families  of  valid  linear  inequalities  (cutting  planes)  have  been  devel- 
oped for  the  two-index  vehicle  flow  formulation  (see  Naddef  & Rinaldi  [25]  for 
a survey).  We  will  be  interested  in  the  following  inequalities: 
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• The  fractional  capacity  (FC)  inequalities: 

x(6+(S))> ^ (SCVc). 

• The  subtour  elimination  (SE)  inequalities: 

x{5+(S))>  1 (S  CVc). 

• The  generalized  large  multistar  (GLM)  inequalities  (see  G 


c(S+{S))  > ^ (*  + X]  q^Xii  + ^J*)) 

^ ies  j&vc\s 


• The  knapsack  large  multistar  (KLM)  inequali  etfchford  et  al.  [21]): 

0)  (SUC),  (9) 

inequality  I]isVc  < P is 


^+(-u)  > nJ2(ai+ 

j£Vc\S 


ies 


where  a > 0 and  /3  > 0 are  such 
valid  for  the  0-1  knapsack  poly 


KP(Q,  q)  : = com  j.y  G {0, 1}”  : ^ q,:y%  < Q j . (10) 

Obviously,  the  RC  ^equalities  dominate  the  FC  and  SE  inequalities,  and 
the  GLM  inequalities  dominate  the  FC  inequalities.  It  is  also  not  difficult  to 
see  that  the  KLM  inequalities  include  the  GLM  and  SE  inequalities  as  special 
cases.  In  general,  there  are  no  other  dominance  relations. 

2. 2^  Single-  and  two- commodity  flow  formulations 

The  jfirst  single-commodity  flow  formulation,  that  we  call  “SCF1”,  was  pre- 
sented by  Gavish  & Graves  [12].  A continuous  variable  ,ftJ  is  defined  for  each 
G A,  representing  the  total  load  (if  any)  carried  along  the  arc  ( i,j ).  One 
then  replaces  constraints  (4)  in  the  two-index  vehicle  flow  formulation  with: 
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(*  e Kc) 

((i,j)eA). 


(11) 

(12) 
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The  constraints  (11)  ensure  that  each  customer  i receives  the  demand  of  qi.  The 
constraints  (12)  are  just  bounds  on  the  / variables. 

Gavish  [13]  proposed  to  strengthen  SCF1  by  replacing  the  bounds  (12)  with 
the  stronger  bounds 

qjXij  < fij  <{Q-  qi)xij  (( i,j ) £ A). 

We  call  this  strengthened  formulation  “SCF2” . 

Gouveia  [15]  used  Hoffman’s  circulation  theorem  [18]  to  proj^  ble 

regions  of  the  LP  relaxations  of  SCF1  and  SCF2  into  the  spaije  of  th^z^variables. 
The  projection  of  SCF1  is  given  by  the  out-degree  equations  (2),  the  in-degree 

if.  "''The 


e projection  of 
ities  (8)  in  place  of 


equations  (3),  the  FC  inequalities  (6)  and  non-negativ: 

SCF2,  as  expected,  is  stronger,  satisfying  the 
the  FC  inequalities. 

In  Baldacci  et  al.  [2] , a two-commodit|'.’fliQw  fbjMiulation  is  presented  for  the 
case  in  which  the  costs  Cij  are  symmetric.  Letcfiford  & Salazar  [22]  show  that, 
in  this  case,  the  LP  relaxation  of  their  formulation  gives  the  same  lower  bound 
as  that  of  SCF2. 

2.3.  Multi- commodity  flow  formulations 

The  first  multi-commodsitf1  flow  formulation,  that  we  call  “MCFla”,  was  pre- 
sented in  Garvin,  et  al.  [11].  A binary  variable  fk  is  defined  for  each  k £ Vc 
and  each  (jf9|\G  A,  taking  the  value  1 if  and  only  if  a vehicle  traverses  (i,  j) 
on  the  wayrrom  the  depot  to  k.  The  formulation  is  then  obtained  by  replacing 
in  the  two-index  vehicle  flow  formulation  with: 


fk(5+(0))  = fk(S~(k))  = l 

(k  £ Pc) 

(13) 

fk(6-(0))  = fk(S+(k))  = 0 

(k  £ Vc) 

(14) 

fk(6+(i))  = fk(5~(i)) 

(k,i  £ Vc  : i k ) 

(15) 

4>° 

(k  £ Vc,  (i,  j)  £ A) 

(16) 

'y  Qkfij  — QXij 

(C hj)£A ). 

(17) 
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The  single-source  multi-commodity  flow  theorem  of  Papernov  [26]  implies  that 
the  LP  relaxations  of  MCFla  and  SCF1  are  of  equal  strength. 

Gavish  [13]  proposed  an  alternative  formulation,  that  we  call  “MCFlb”.  It  is 
obtained  by  replacing  (17)  with  the  following  constraints: 


£ qkfk(S+(i))<Q~q. , 
kevc\{i} 

fk 


(*  e Vc) 


fij  < Xij  (fc  e Vc,  (i,j)  £ A).  (1 

amts  £ 


A 

onstran 


It  follows  from  the  max-flow/min-cut  theorem  that,  if  the  constraints  flS J are 
dropped  from  MCFlb,  then  the  projection  into  x-space  is  given  by  (2).  (3),  (7) 
and  non-negativity.  No  similar  projection  result  is  known  for  MCFlb  itself. 

Letchford  & Salazar  [22]  presented  a different  MCF  formulation,  with  two 
commodities  per  customer.  For  each  arc  (i,j)  and  customer  fc,  the  variable  /*) 
is  defined  as  before,  but  there  is  now  also  a binary  variable  gfj , taking  the  value 
1 if  and  only  if  a vehicle  traverses  (*,  j)  on  the  way  from  k to  the  depot.  We 
then  replace  the  constraints  (4)  in  the  two-index  vehicle  flow  formulation  with 
the  constraints  (13)— (16),  together  with: 


(urvon  o))  = i 

4/(J+(0))  = 0 
(S+(i))  = gk(S~(i)) 
9%  > 0 
fij  + 9ij  < 
{ftS+m  + g^S-ik)))  <Q-qi 


(fc  £ Vc)  (20) 

(fc  e Vc)  (21) 

(k,i  £ Vc  : i ^ k)  (22) 

(fc  £ Vc,  (i,  j)  £ A)  (23) 

[k  £ Vc,  (i,  j)  £ A)  (24) 

( i £ Ve).  (25) 


call  this  formulation  “MCF2a”.  Note  that  the  depot  is  either  the  source 
the  sink  of  every  commodity.  The  above-mentioned  result  by  Papernov  [26] 
then  implies  that,  if  the  constraints  (25)  are  dropped  from  MCF2a,  then  the 
projection  into  x-space  is  again  given  by  (2),  (3),  (7)  and  non-negativity.  No 
similar  projection  result  is  known  for  MCF2a  itself. 
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2.f.  Set  partitioning  formulations 

We  will  also  need  the  following  set  partitioning  (SP)  formulation,  due  to 
Balinski  & Quandt  [4].  Let  12  denote  the  set  of  possible  routes  for  a single 
vehicle,  and  let  zr  for  each  r G 12  be  a binary  variable  taking  the  value  1 if  and 
only  if  that  route  is  used.  Define  the  constant  a*r  for  each  customer  i and  jafeitel 
r,  taking  the  value  1 if  i is  served  by  r,  and  0 otherwise.  Finally  let  cf 
the  cost  of  route  r.  Then  the  SP  formulation  is: 

min  crzr 

r(E  f2 

s.t.  airzr  = 1 

zr  G {0, 1} 

Since  the  number  of  variables  in  this  forffi  fctfoftAan  be  exponential  in  n, 
column  generation  is  necessary.  Unfortunately,  tHJfpricing  subproblem  is  easily 
shown  to  be  strongly  A/’'P-hard.  Agarwal  et  al.  [1]  solve  it  via  integer  program- 
ming. Foster  & Ryan  [9]  noted  that  pricing  becomes  easier  if  one  enlarges  the  set 
17  by  allowing  routes  in  which  the  vehicle  is  permitted  to  visit  customers  more 
than  once  (now  called  non- elementary  routes).  Pricing  can  then  be  performed 
in  pseudo-polynomial  time^by  dynamic  programming.  See,  e.g.,  Martinelli  et 
al.  [24]  for  details 

Letchford  and  Salazar  [22]  prove  the  following: 

/N  7 

1 elementary  routes  are  used,  the  projection  of  the  LP  relaxation  into 
ace  satisfies  all  KLM  inequalities. 


ain,  when  elementary  routes  are  used,  the  LP  relaxation  is  at  least  as 
strong  as  those  of  all  of  the  SCF  and  MCF  formulations  mentioned  in  the 
previous  two  subsections. 

• If,  however,  non-elementary  routes  are  permitted,  then  the  only  KLM 
inequalities  which  are  satisfied  by  the  projection  are  those  in  which  ay  < (3 
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is  valid  for  the  general  integer  knapsack  polytope 

conv  {y  G Z™  : E QiVi  < Qj- 
ievc 

These  less  general  KLM  inequalities  still  include  the  GLM  inequalities  a: 
a special  case,  but  no  longer  include  the  SE  inequalities. 


3.  Stronger  Multicommodity  Flow  Formulations 


In  this  section,  we  present  four  new  multi-commodity  flow  formulations 
each  of  which  satisfies  the  SE  inequalities  (7)  and  GLM.fiAequalities  (8).  The 
ones  presented  in  Subsection  3.1  dominate  SCF1,  SCFSOlCFla  and  MCFlb.  The 
one  presented  in  Subsection  3.2  also  dominates  MCF2a.  The  one  presented  in 
Subsection  3.3  is  designed  especially  for  instances  with  symmetric  costs. 

3.1.  Strengthening  MCFla  and  MCFlb 

In  this  subsection  we  will  need  the,  foll&fin^lemma. 

Lemma  1.  The  LP  relaxation  of  formulation  MCFlb  satisfies  the  equations 

j (j  eVc,ieV\{j}).  (26) 


Proof.  Let  j G Vc  be  fixed.  From  (3)  we  have  x(5~(j))  = 1,  from  (13)  we  have 
/J(<5-(j))  = 1,  and  from  (19)  we  have  f^  < Xij  for  all*  G V \ (j).  The  only 
way  for  these  to  all  hold  simultaneously  is  for  (26)  to  hold.  □ 

The  flowing  proposition  introduces  a class  of  valid  inequalities. 

Proposition  1.  All  (integer)  solutions  to  formulation  MCFlb  satisfy  the  follow- 
ing inequalities: 


E Qkfij  <{Q~  Qi > 

k£Vc\{i} 


1 3 


{{hi)  e A). 


(27) 


Proof.  If  the  vehicle  traverses  the  arc  ( i , j ),  then  it  must  have  already  delivered 
a demand  of  qi  to  customer  i.  □ 
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Our  first  new  formulation,  which  we  call  MCFlc,  is  obtained  from  MCFlb  by 
replacing  the  constraints  (18)  with  inequalities  (27).  The  following  two  propo- 
sitions state  that  MCFlc  has  some  desirable  properties. 

Proposition  2.  The  LP  relaxation  of  MCFlc  satisfies  the  SE  inequalities 
and  the  GLM  inequalities  (8). 

Proof.  As  mentioned  in  Subsection  2.3,  the  max-flow  / min-' 
implies  that  the  SE  inequalities  are  satisfied.  Now,  for  a give: 
given  commodity  fe  G S',  the  flow  equations  (13)— (15)  imply  t 

ck  ( £— 


f\5~{S))  = fk{5+{S)) 


Multiplying  these  equations  by  qk  and  summin, 


G S,  we  obtain 


E«ufc(<nu)  = E^/fc(^<  ))+^(s). 

fees  fceS  \ 

<?—  (S')  iir 


Now,  the  constraints  (27)  for  all  ) imply  that  the  left-hand  side  of 

(28)  is  no  larger  than 

'y  \ ( Q — qi)^ij- 

On  the  other  hand.  tetreunts  (26)  for  all  (*,  j)  G 5+(S),  together  with 

non-negativity,  iia#  ^Ihe  right-hand  side  of  (28)  is  no  smaller  than 

E QjXij  + q{S). 

(i,j)es+(s) 


educe  that  the  relaxation  of  MCFld  satisfies 

E {Q-Qi)xij>  E <Hxij+q(S)’ 

ich  is  equivalent  to  the  GLM  inequality  (8)  for  the  given  S. 


SP 

Y'  (/ 

-ti 

</ 

X3-n 

QJ-n 

o>m 
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Proposition  3.  The  LP  relaxation  of  MCFlc  is  stronger  than  that  of  SCF1,  SCF2 
and  MCFla,  and  at  least  as  strong  as  that  of  MCFlb. 
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Proof.  The  LP  relaxation  of  MCFlc  satisfies  the  constraints  (2),  (3)  (7)  and 
(8).  The  fact  that  it  is  stronger  than  the  LP  relaxations  of  SCF1  and  SCF2  then 
follow  from  the  result  of  Gouveia  [15]  mentioned  in  Subsection  2.2.  It  is  also 
stronger  than  the  LP  relaxation  of  MCFla  since,  as  mentioned  in  Subsection  2.3, 
that  relaxation  is  identical  to  the  one  of  SCF1.  To  show  that  the  LP  relaxg|M 
of  MCFlc  is  at  least  as  strong  as  that  of  MCFlb,  it  suffices  to  show  that  it  i 
the  inequalities  (18).  To  this  end,  let  * G Vc  be  fixed.  Sum  the  inec 
over  all  arcs  entering  i to  obtain 

E <Q  - 53  qkXk 

kev c fceVc\{i} 


Together  with  the  equation  (13)  for  i = k,  this  imp 

55  Qkfk(S-{i))<Q-qi J 

fcCVc\T} 

The  equations  (15)  then  imply 


E f(S+(i )) 

fceVcUd 

This  dominates  the  inequality 


E 

fcevc\{i} 
the  given  i. 


Qk-^ki  • 


ineq 

tidn. 


Our  computational  results  (Section  5)  show  that,  in  fact,  the  LP  relaxation 
of  MCFlc  is  stronger  thaB^that  of  MCFlb,  and  is  also  stronger  than  the  relaxation 
in  a;-space  defined  by  the  out-degree  equations  (2),  the  in-degree  equations  (3), 
the  SE  inequalities  (7,),  the  GLM  inequalities  (8),  and  non- negativity. 

Now,  in  the  proof  of  Proposition  3,  we  showed  that  the  inequalities  (18) 
are  rMujraairtTbeing  implied  by  the  other  constraints  in  MCFlc.  Interestingly, 
however^Mey  can  be  strengthened  (lifted)  to  obtain  a non-redundant  family  of 


inequalities.  These  stronger  inequalities  are  presented  in  the  following  proposi- 
tic 


Proposition  4.  All  (integer)  solutions  to  formulation  MCFlc  satisfy  the  follow- 
ing inequalities: 

E Qk(f(S+m  + fk(6+(i)))  <Q-qi  (ie  Vc).  (29) 
fcevc\{i} 
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Proof.  Consider  the  vehicle  that  services  customer  i.  The  two  terms  on  the 
left-hand  side  represent  the  total  demand  delivered  by  the  vehicle  before  arriving 
at  customer  i and  after  leaving  customer  i,  respectively.  □ 


Adding  the  inequalities  (29)  to  MCFlc,  we  obtain  a formulation  that  we 
call  “MCFld”.  (For  ease  of  reference,  we  present  MCFld  in  its  entiretyArSthe 
Appendix.)  The  computational  results  in  Section  5 show  that  the  LP  refaxation 
of  MCFld  is  stronger  than  that  of  MCFlc,  which  shows  that  the  ineqiialitie 
are  not  implied  by  the  other  linear  constraints  in  MCFld. 

3.2.  Strengthening  MCF2a 

Now  we  turn  our  attention  to  MCF2a.  We  will  sfcqw  tha/  MCF2a  can  be 
strengthened  to  obtain  a formulation  that  also  dominates  MCFld,  the  strongest 
of  the  formulations  given  in  the  previous  subsnctioN  Our  starting  point  is  the 
following  proposition: 

Proposition  5.  All  (integer)  solutiq^s  to^rrr^ulation  MCF2a  satisfy  the  follow- 
ing constraints: 

(k,i  G Vc  : i k).  (30) 

-hand  sides  are  each  equal  to  1 if  customer  i is 
same  route,  and  0 otherwise.  □ 

serve  that,  using  Lemma  1,  constraints  (27)  can  be  written 


Proof.  The  left-ha: 
served  before 

The  next  stq 
as 


Qkftj  <{Q~  q:  ~ qj)a 


(ihj)  e A), 


kevc\{i,j} 


following  proposition  shows  that  these  inequalities  can  then  be  strength- 
ened, using  the  ^-variables. 

Proposition  6.  All  (integer)  solutions  to  formulation  MCF2a  satisfy  the  follow- 
ing constraints: 

Y (fij  + 9ij)  <(Q-qi~  qj)xij  (( i,j ) G A).  (31) 

kevc\{  i,j} 
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Proof.  If  a vehicle  traverses  the  arc  ( i , j ),  then  it  must  have  already  delivered  a 
demand  of  qi  to  customer  i,  and  be  about  to  deliver  qj  to  customer  j.  Therefore 
the  total  load  of  the  vehicle  dedicated  to  the  other  customers  on  the  route  cannot 
exceed  Q — qi  — qj. 

We  now  show  that  the  constraints  (25)  can  be  discarded  when  (30)  an 
are  considered.  To  this  end,  we  need  the  following  lemma: 

Lemma  2.  The  LP  relaxation  of  formulation  MCF2a  satisfies 


fii  = Uh  = xk 


ihj  eVc:i^  j) 


Proof.  Similar  to  the  proof  of  Lemma  1. 

Proposition  7.  If  the  constraints  (30)  and  (31 
constraints  (25)  become  redundant. 


Proof.  Use  (21)  and  (32)  to  write 

qpXij 

nd  use  (2)  to  obtain: 


H 9k  If!)  + g{ 
feev0\{*} 


Sum  these  over  all  j 


((*>j)  e A- 


(i))+gk(S+(f)))<Q-qi  (i  e Vc). 


Reverse  the  roles  of  i and  k.  and  use  (15)  and  (30)  to  obtain  (25). 


□ 


the 


Aetorai'nMy,  we  add  constraints  (30)  and  (31)  to  MCF2a,  and  delete  the 
jfjedundant  constraints  (25),  to  obtain  what  we  call  “MCF2b”.  (We  present  MCF2b 
in  its  entirety  in  the  Appendix.)  The  following  proposition  shows  that  MCF2b  is 
the  strongest  of  all  the  MCF  formulations  considered  so  far. 


Proposition  8.  Let  (x*,f*,g*)  be  a solution  to  the  LP  relaxation  of  MCF2b. 
Then  (x* , f*)  is  a solution  to  the  LP  relaxation  of  MCF  Id. 
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Proof.  The  feasible  region  of  the  LP  relaxation  of  MCFld  is  defined  by  the 
constraints  (2),  (3),  (13)— (16),  (19),  (27)  and  (29).  Constraints  (2),  (3)  and 
(13) — (16)  are  already  present  in  MCF2b.  Constraints  (19)  are  dominated  by 
(24).  Constraints  (27)  are  dominated  by  (31).  Finally,  the  fact  that  constraints 
(29)  are  satisfied  by  (x *,/*)  follows  from  Proposition  7 together  with  theyract; 
that  MCF2b  contains  the  equations  (30). 

Corollary  1.  The  LP  relaxation  of  formulation  MCF2b  satisfies  th4^Ejjjgi 
ities  (7)  and  the  GLM  inequalities  (8). 

Our  computational  results  (Section  5)  indicate  that  t^^^syer  bound  given 
by  MCF2b  is  signficantly  stronger  than  the  lower  bounds  given  bymither  MCFld  or 
MCF2a.  It  is  therefore  the  strongest  known  formuljjfiol|^)ased  on  commodity-flow 
variables. 

3.3.  An  alternative  formulation  with  fewer  vanaOfes 

To  finish  this  section,  we  present  a ‘hybrid’  formulation,  which  attempts  to 
‘mimic’  formulation  MCF2b,  but  use  only  one  commodity  per  customer,  instead 
of  two.  We  will  show  that,  when  a CVRP  instance  has  symmetric  costs,  the 
lower  bound  from  this  hybrid  formulation  is  just  as  good  as  the  one  from  MCF2b. 

The  basic  idea  of  the  ‘IWMjia’  formulation  is  to  take  the  formulation  MCF2b 
(given  in  the  Appendix) , and  simplify  it  by  using  aggregated  flow  variables  fk-  = 
fij  +&  f°r  * U and  (i,  j)  £ A.  (Note  that  these  aggregated  variables  have 
no  natural  physical  interpretation.)  The  constraints  in  the  resulting  formulation 
follows: 


/fe(<5+(0))te/fe(<r(fc))  =2 
fk{6~m  = fk(S+(k))=0 
fk(6+(i))mfk(6-(i)) 

fii>  o 


fij  + fji  < Xij  + x 


Jl 


Y 9k  ( ftj  + fji)  < (Q  - Qi  ~ Qj ) (xij  + Xji ) 

keVa\{ij} 


(k  e vc ) 

(k  £ Vc ) 

(k,  i £ Vc  : i 7^  k) 

( i,j,k£  Vc  :i/j) 
(i,j,  k £ Vc  : i < j) 
(' i,j  £Vc:i<j) 
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fk(5+(i))  = f(S+(k))  (k,ieVc:i^  k). 

We  call  this  formulation  “MCF3”.  Although  MCF3  is  weaker  than  MCF2b  by 
construction,  it  does  have  some  interesting  properties.  These  are  given  in  the 
following  theorem  and  corollaries. 


Theorem  1.  If  the  costs  are  symmetric,  i.e.  Cij  = Cji  for  all  ( i,j ) E A, 

the  lower  bounds  obtained  by  solving  the  LP  relaxations  of  MCF2b  and  MCF3'&j;e 

A 

equal. 

(TV 

Proof.  Let  (x,  f)  be  any  feasible  solution  to  the  LP  relaxation  of  MCF3.  We 
can  construct  a feasible  solution  to  the  LP  relaxation  of  MCF2b,  say  ( x,f,g ), 
using  the  following  mapping: 


Xij  — (,X{j  T Xjfjf 2 

fii = = 


m 


A,  k E Vc). 


Moreover,  if  the  costs  are  symmetri4^,hen  the  cost  of  ( x , /,  g ) is  the  same  as 
that  of  ( x,f ).  This  shows  that,  when  costs  are  symmetric,  the  lower  bound 
from  MCF2b  cannot  be  better  than  the  one  from  MCF3.  But  it  cannot  be  worse, 
since  MCF3  is  weaker  than  MCF2b  by  construction.  □ 

Corollary  2.  The  LP  relaxation  of  MCF3  satisfies  the  SE  inequalities  (7)  and 
the  GLM  inequalities  (8). 

Proof.  The  in-  and  out-degree  equations  imply  that  aj(<S+(‘S'))  = x(S~(S))  for 
all  S / Vv  It  fellows  that  the  SE  inequalities  can  be  written  as 

a:(<S+(S))-l-a:(«r(S))>2  (SCVc), 

e GLM  inequalities  can  be  written  as 

(<U(A))  + x(5~(S))  > ^ (&  + qi(Xii  + xa))  ( S £ vc)- 

^ ies  j£Vc\S 

In  this  form,  the  inequalities  are  symmetric  (that  is,  for  all  (i,j)  E A,  the 
coefficients  of  Xij  and  Xji  are  equal).  The  result  then  follows  from  Corollary  1 
and  Theorem  1.  □ 
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Corollary  3.  The  LP  relaxation  of  MCF3  is  stronger  than  those  of  SCF1,  SCF2 
and  MCF1  a. 


^property  of 


Our  computational  results,  given  in  Section  5,  suggest  that  the  LP  relaxation 
of  MCF3  may  be  stronger  than  those  of  MCFlb,  MCFlc  and  MCF2a  as  well.  On  tl 
other  hand,  they  show  that  there  is  no  dominance  relationship  between 
relaxations  of  MCF3  and  MCFld. 

4.  Stronger  Bounds  in  Pseudo-Polynomial  Time 

In  this  section  we  show  that,  if  we  are  willing  to 
having  only  a polynomial  number  of  variables  and  constraints,  then  we  can 
strengthen  both  MCFld  and  MCF2b  even  further.  In  Subsection  4.1,  we  strengthen 
them  by  adding  certain  inequalities  derived  from  fac^s  of  the  0-1  knapsack  poly- 
tope, and  show  that  the  continuous  relaxations  of  the  resulting  formulations 
satisfy  all  KLM  inequalities.  In  Subsection  4.2,  we  present  alternative  strength- 
ened formulations,  of  the  same  quality,  which  have  additional  variables,  rather 
than  constraints.  We  then  show  thatfflae^ontinuous  relaxations  of  these  latter 
formulations  can  be  solved  in  pseudd^polynomial  time,  via  column  generation. 

4-1.  Valid  inequalitiestfroin  the  knapsack  polytope 

The  followingjigmma  introduces  exponentially-large  families  of  valid  inequal- 
ities that,  in  theory,  could  be  used  to  strengthen  further  the  formulations  MCFld 
and  MCF2b. 

Lemma  3.  Let  aTy  < f3  be  any  valid  inequality  for  the  0-1  knapsack  polytope 
(10),  with  a > 0 and  (3  > 0.  Then  all  (integer)  solutions  to  formulation  MCFld 
isfy  the  inequalities 


^2  akfij  < (P-  ati)xij 


(0J)  e A), 


fceVcUd 

and  all  (integer)  solutions  to  formulation  MCF2b  satisfy  the  inequalities: 
ak  Uij  + 9ij ) < (/?  - on  - aj)xij  (0 i,j ) £ A). 

kev0\{i,j} 


(33) 


(34) 
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Proof.  Let  ( x,f ) be  a feasible  integer  solution  to  MCFld.  If  aty-  = 0,  the 
inequality  (33)  holds  triviality.  So,  suppose  that  Xy  = 1,  i.e.,  a vehicle  traverses 
the  arc  (*,  j).  Now  consider  the  vector  y e (0,  l}ra  obtained  by  setting  yk  to  1 
if  and  only  if  commodity  k is  on  the  vehicle  just  before  the  vehicle  arrives  ajj 
node  i.  Then  y must  be  an  extreme  point  of  the  polytope  (10),  and  thereto 
satisfies  aTy  < (3.  Moreover,  jh  must  equal  1,  which  implies 

^2  akyk  < (p  - 

k£Vc\{i} 


,tisfies  the 


Now,  for  k G Vc  \ {i},  we  must  have  yk  = fij-  This  impli^ 
inequality  (33). 

The  proof  for  MCF2b  is  similar.  The  only  differences  are  that  (i)  we  set  yk  to 
1 if  and  only  if  commodity  k is  on  the  vehicle  when  it  leaves  the  depot  and  (ii) 
for  k e Vc  \ {i,  j},  we  must  have  yk  = fij  + gfp  □ 

We  have  the  following  result: 


Theorem  2.  If  all  non-redundant  inequalities  of  the  form  (33)  are  added  to 
MCF1  d, , the  continuous  relaq^Upn  of  vjtfe  resulting  formulation  satisfies  the  KLM 
inequalities  (9).  The  sdme  holds  if  all  non-redundant  inequalities  of  the  form 
(34)  are  added  to  MCFZb 

Proof.  To  show  that  the  relaxation  of  MCFld  with  (33)  added  satisfies  the 
KLM  ineqhahpes  (9),  just  use  the  proof  of  Proposition  2 with  <Xi  and  (3  taking 
the  pbfce'of  ty  and  Q.  and  (33)  taking  the  place  of  (27).  To  show  the  same  result 
for  MCF2b  with  (34)  added,  just  note  that  the  inequalities  (34)  are  stronger  than 
inequalities  (33).  □ 

4- f.  Alternative  formulations  solvable  by  column  generation 

The  models  in  the  previous  subsection  have  a drawback  from  a computa- 
tional point  of  view.  Not  only  are  the  constraints  (33)  and  (34)  exponential  in 
number,  but  their  coefficients  are  given  only  implicitly.  To  overcome  this  draw- 
back, we  present  alternative  formulations  that  are  more  explicit,  but  yield  lower 
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bounds  of  the  same  quality.  These  formulations  have  an  exponential  number  of 
variables  rather  than  constraints. 

Let  P be  the  set  of  all  possible  loading  patterns  of  a single  vehicle.  That 
is,  each  member  of  P is  a subset  of  Vc  whose  total  demand  does  not  exceed  Q. 
For  each  arc  (i,  j)  and  for  each  p G P,  let  X1’,-  be  a binary  variable,  taking 
value  1 if  and  only  if  a vehicle  arrives  at  node  i carrying  loading  pattern 
then  traverses  arc  (i,  j).  (This  implies  of  course  that  i and  j are  ir 
the  model  MCFld  is  enlarged  by  adding  the  constraints 


4 


4 


e A:  i,j  g Vc) 
((i,  j)  € A:  i,j  G 14,  fc 


Xij  ~ E^ 

peP:  i,jep 

/«  = E 

pEP:  i,j,k£p 

Now,  for  each  arc  (i,  j)  and  for  each  p G P,  let 
the  value  1 if  and  only  if  a vehicle  that  depar 
pattern  p goes  on  to  traverse  the  arc 
in  p).  Then  the  model  MCF2b  is  enlai 


XH  = 

p£P:  ijtzp 

4+4=  E M 

p&P:  tj,k£p 


inary  variable,  taking 
rom  the  depot  with  loading 
mplies  again  that  i and  j are 
dding  the  constraints 

(37) 

(38) 


G A:  i,j  G Vc) 

(( i,j ) e A'-  i,j  eVc,ke  Vc\{i,j}). 


Observe  that  the  members  of  P correspond  to  extreme  points  of  the  0- 
1 knapsack  polytope  (10).  Accordingly,  we  call  the  strengthened  formulations 
“MCF1K”  and  ‘:MCF2# , respectively.  The  following  two  theorems  show  that  these 
formulations  have  two  desirable  properties. 

Theorem  3.  The  LP  relaxations  of  formulations  MCF1K  and  MCF2K  satisfy  the 
KLM  inequalities  (9). 

Proof.  Let  a and  /3  be  as  in  Lemma  3.  From  the  definition  of  P it  follows 
that  ak  — @ f°r  a+  P £ P-  Therefore,  for  all  t G kc  and  all  p G P such 

that  i G p,  we  have  Y^kep\{i}  ak  < /3  — ai.  This  implies 


E E °*  Uf<G8-««)  E 


A?.- 
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for  all  (i,j)  G A.  Rearranging  the  left-hand  side,  and  using  (35)  to  simplify  the 
right-hand  side,  we  get 


^2  Q-k  ^2  A <{/3-  a.i)xij  (( i,j ) G A:  ij  G Vc). 

fcgVc\{j}  pGP:  ij'jfcGp 

These  constraints  together  with  (36)  then  imply  (33).  In  a similar  way, 
straints  (37)  and  (38)  imply  (34).  The  result  then  follows  from  Theorem  2. 
□ 

Theorem  4.  The  LP  relaxations  of  MCF1K  and  MCF2K  can  be  pseudo- 

polynomial time. 

Proof.  Suppose  we  have  solved  a restricted  master  pr  em  associated  with 
MCF1K,  i.e. , an  LP  obtained  from  MCF1K  by  relaxing  the  integrality  condition  and 
replacing  P with  a small  (polynomial-sized)  subset  P'  C P.  Let  pij  and  7r£)  be 
the  optimal  dual  prices  for  (35)  and  (36),  respectively.  For  a given  ( i,j ) G A, 
there  exists  a column  A with  negati;  ^£ed%zed  cost  if  and  only  if 


max  • 


yk<Q-qi-  qj , y g {0,  i}n  2 [•  > ptJ. 

k£V0\{i,j} 

This  is  a 0-1  knapsack  problem,  that  can  be  solved  in  0(nQ ) time  via  dynamic 
programming  [t>L  iiie  pricing  problem  for  MCF2K  is  similar.  Now  note  that  pric- 
ing in  an  LP  is  equivalent  to  separation  in  the  dual  of  the  LP.  The  desired  results 
then  follow  from  theq>olynomial  equivalence  of  separation  and  optimization  [16]. 
□ 


We  rejliark  that  this  is  the  first  time  that  an  LP  relaxation  of  the  CVRP 
fajt&pbn  found  that  can  be  solved  in  pseudo-polynomial  time,  yet  satisfies  all 
he  KLM  inequalities.  (As  mentioned  in  Subsection  2.4,  the  LP  relaxation 
of  the  SP  formulation  with  elementary  routes  satisfies  the  KLM  inequalities, 
but  is  strongly  TVP-hard  to  solve.)  Moreover,  no  pseudo-polynomial  separation 
algorithm  is  known  for  the  KLM  inequalities  themselves. 
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f.3.  Comparison  with  set  partitioning 

The  strongest  of  our  MCF  formulations  is  MCF2K.  A natural  question  is  how 
the  lower  bound  associated  with  MCF2K  compares  with  the  lower  bound  asso- 
ciated with  the  SP  formulation  given  in  Subsection  2.4.  The  answer  to  this 
question  depends  on  whether  or  not  the  set  fi  is  permitted  to  contain 
that  correspond  to  noil-elementary  routes.  As  mentioned  at  the  end 
section  2.4,  the  LP  relaxation  with  non-elementary  routes  permitted 
satisfy  all  KLM  inequalities.  Thus,  the  associated  bound  does  not  dominate  the 
bound  from  MCF2K.  The  following  theorem  settles  the  que^ion  for  the  case  m 
which  only  elementary  routes  are  permitted: 
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Theorem  5.  If  z*  is  a solution  to  the  LP  relaxat 
then  there  exists  a solution  (x*,  /*,  <?*,  fi*)  to 
has  the  same  cost. 


of  the  SP  formulation, 
relaxation  of  MCF2K  that 


Proof.  For  each  arc  (i,  j)  G A and  each  route  r G fi,  let  bijr  be  a binary  con- 
stant which  takes  the  value  1 if  and  only  if  route  r traverses  arc  (i,  j).  Also,  for 
any  {i.j)  G A,  any  r G fi  and  any  ch^omer  node  k G Vc,  let  dijkr  (respectively, 

d'nkr)  be  a binary  constant  which  takes  the  value  1 if  and  only  if.  in  route  r, 

/\  / 

/ v v 

• vertex  k is  visited  and 


• the  arc  w,  ff  is  traversed  on  the  way  from  the  depot  to  k (respectively,  on 
the  way  from  k to  the  depot). 

Finally,  for  any  (i,j)  G A,  r G fi  and  loading  pattern  p G P,  let  Ujrp  be  a 

binary  constant  which  takes  the  value  1 if  and  only  if  a vehicle  following  route 

r departs  from  the  depot  with  loading  pattern  p and  then  goes  on  to  traverse 
ay  (i.  j).  The  desired  quadruple  (a "* , f* , g* , p,*)  is  then  created  by  setting: 

xij  t0  Ere  fi  KrK  for  a11  (hi)  G A 
/*/  to  Eren  dijkrZ*  for  all  (i,j)  G A and  all  k G Vc; 

9ij  *0  dijkrZr  for  a11  (hi)  G A alld  a11  k G VA 
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• MiJ  to  J2ren  tijrpZ*  for  a11  (hj)  G A and  P G P- 

This  quadruple  has  the  same  cost  as  z* , from  the  definition  of  cr  in  the  objective 
of  the  SP  formulation.  One  can  check  that  it  also  satisfies  all  of  the  linear 
constraints  in  the  formulation  MCF2K,  i.e.,  the  constraints  (2),  (3),  (13) — (16) 
(20)-(24),  (31),  (30),  (37)  and  (38). 

Theorem  5 implies  that  the  lower  bound  from  MCF2K  is  dominated 
bound  from  the  SP  formulation  with  elementary  routes.  We  stre  lat 

the  former  bound  can  be  computed  in  pseudo-polynomial  thr  efeas  the 
latter  bound  cannot  (unless  V = MV). 

5.  Computational  Experiments 

In  this  section,  we  report  on  some  computational  experiments.  We  stress 
from  the  outset  that  the  goal  of  these  experiments  was  not  to  solve  large-scale 
CVRP  instances  to  proven  optimality,  but  ratDer  to  establish  dominance  re- 
lations between  the  lower  bounds  obtained  when  solving  the  LP  relaxation  of 
various  formulations.  (The  development  of  a viable  exact  algorithm  for  the 
CVRP  based  on  formulation  MCF2K  may  be  the  topic  of  a future  paper.)  We 
found  that,  for  this  purpose,  it  was  sufficient  to  use  small  instances  with  n = 16. 
The  advantage  of  using  these  instances  is  that  the  RC  inequalities  (4)  can  be 
enumerated  and  added  to  the  LP  relaxation  if  desired.  (No  efficient  separation 
procedure  is  known  for  the  RC  inequalities.) 

We  created  both  asymmetric  and  symmetric  instances.  In  the  asymmetric 
instanc%,.  the  costs  Cij  were  randomly  generated  in  [0,500].  In  the  symmet- 


ric instances,  the  costs  were  obtained  by  computing  the  Euclidean  distance 
between  locations  randomly  distributed  in  the  square  [0,500]  x [0,500],  except 
the  depot,  which  was  located  in  the  center  of  the  square.  We  created  instances 
with  general  demands  (random  integers  in  the  range  [25,33])  and  instances 
with  only  unit  demands.  For  the  instances  with  general  demands,  we  consid- 
ered Q G {100,150,200}.  For  the  instances  with  unit  demands,  we  considered 
Q G {4,6,8}.  This  led  to  twelve  families  of  instances,  and  for  each  family  we 
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Table  1:  Average  ratios  for  relaxations  in  .x'-spa.ce. 


Type 

SE 

FC 

GLM 

RC 

FC+SE 

GLM+SE 

GLM+RC 

A-G-100 

42.47 

74.76 

80.09 

91.61 

74.96 

80.10 

91.82 

A-G-150 

64.10 

82.56 

84.57 

89.88 

83.03 

84.77 

89.88 

A-G-200 

72.64 

81.35 

82.14 

96.62 

82.17 

82.70 

96.62 

A-U-4 

54.37 

83.99 

87.92 

94.30 

84.33 

87.92 

94.35 

A-U-6 

72.12 

85.85 

87.34 

97.65 

86.62 

87.79 

97.65 

A-U-8 

87.37 

93.76 

94.33 

98.05 

94.98 

95.25 

98.05V 

S-G-100  61.53 

S-G-150  76.83 

S-G-200  82.31 


81.52 

82.85 

79.70 


87.48 

86.91 

82.07 


98.70  85.42 

99.62  90.95 

99.64  90.25 


89.65 

92.95 

91.38i 


S-U-4  70.01 

S-U-6  81.17 

S-U-8  89.51 


85.20 

82.43 

82.84 


90.76 

85.43 

84.71 


99.98 

99.95 

100.00 


90.66 

91.94 

95.44 


generated  20  instances.  The  instance  generator  i the  formulations  were  im- 
plemented in  FICO  Xpress  Mosel  3,  and  the  source  code  is  available  to  readers 
on  request  to  the  authors. 

Table  1 gives  the  results  for  seven  LP  relaxations  that  only  involve  x- 
variables.  The  first  column  descries  the  instance  type,  and  the  remaining 
columns  summarise  the  results  that  we  obtained  when  adding  various  combina- 
tions of  the  SE,  FC,  Gl/M  and  RC  inequalities  to  the  LP  relaxation  consisting 
of  (l)-(3)  and  non- negativity.  Each  figure  is  the  average,  over  20  instances,  of 
the  ratio  between  the  lower  bound  and  the  optimum,  expressed  as  a percentage. 
Table  2^gives  analogous  results  for  six  LP  relaxations  that  involve  /-variables, 
Table  3 does  the  same  for  five  relaxations  that  involve  /-  and  g- variables,  and 
Table  4 does  the  same  for  four  SP  relaxations,  that  involve  x and  2 variables. 
Columns  nSP  and  eSP  refer  to  the  SP  relaxation  with  non-elementary  and  el- 
ementary routes,  respectively.  Columns  nSP+RC  and  eSP+RC  refer  to  the  same 
relaxations  strengthened  with  RC  inequalities  (4). 

The  main  conclusion  from  the  results  in  Table  1 is  that  the  RC  inequalities 
are  the  most  important  inequalities  by  far.  Comparing  Tables  1 and  2,  we  see 
that,  as  expected,  MCFla  gives  the  same  bound  as  the  FC  inequalities.  We 
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Table  2:  Average  ratios  for  relaxations  in  (x , /)-space. 
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3i 


Table  4:  Average  ratios  for  relaxations  in  ( x , z)-space. 


Type 


nSP  nSP+RC  eSP  eSP+RC 


A-G-100  98.56 

A-G-150  93.40 

A-G-200  89.38 


98.72 

93.76 

97.45 


98.59 

94.09 

90.94 


98.72 

94.25 

97.67 


A-U-4  94.76 

A-U-6  92.11 

A-U-8  96.97 


96.86 

98.19 

98.28 


94.97 

93.39 

98.43 


96.92 

98.40 

98.86 


S-G-100  95.05 

S-G-150  90.78 

S-G-200  85.37 


99.39 

99.67 

99.70 


98.46 

98.69 

97.16 


99.76 

99.88 

99.83 


S-U-4  91.91 

S-U-6  86.66 

S-U-8  85.74 


100.00 

99.95 

100.00 


98.29 

97.27 

99.50 


also  see  that  the  four  new  formulations  give  better  bounds  than  MCFla  and 
MCFlb,  and  also  better  bounds  than  the  iSE  and  GLM  inequalities  combined. 
Note  also  that  there  is  no  dominance  between  MCFla  and  MCFlb,  but  MCFlb 
tends  to  perform  poorly.  Moreover,  as  expected,  MCFlc,  MCFld  and  MCF1K  are 
of  increasing  strength.  As  for  MCF3,  jye  see  that  it  appears  to  dominate  MCFlc, 
but  does  not  dominate  MCFld.  Interestingly,  MCF1K  does  not  dominate  MCF3, 
despite  the  fact  that  it  satisfies  all  KLM  inequalities  and  is  weakly  A/”'P-hard  to 
compute.  Also,  MCF1K  gives  the  same  bound  as  MCFld  in  the  unit  demand  case. 

Turning  our  attention  to  Table  3,  we  see  that  MCF2a  consistently  gives  worse 
bounds  than  the  SE  and  GLM  inequalities  combined.  As  expected,  MCF2b  dom- 
inates MCFld  and  MCF2a,  and  gives  the  same  bound  as  MCF3  on  the  symmetric 
instances}  Also  as  expected,  MCF2K  dominates  MCF1K  and  MCF2b.  In  fact,  it  is 


✓^•significa 
the  sami 


significantly  stronger  than  MCF1K  in  all  cases.  On  the  other  hand,  MCF2K  gives 
thfe  same  bound  as  MCF2b  in  the  unit  demand  case.  Note  also  that  using  MCF2b 
or  MCF2K  in  combination  with  the  RC  inequalities  gives  better  results  than  using 
RC  and  GLM  inequalities  in  combination.  The  difference  is  noticeable  especially 
for  the  asymmetric  instances. 

Now  consider  Table  4.  We  see  that,  as  expected,  the  SP  relaxation  with 
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non-elementary  routes  dominates  the  GLM  relaxation.  On  the  other  hand,  it 
is  weaker  even  than  MCFlc  in  some  cases,  despite  being  weakly  A/’T’-hard  to 
compute.  Moreover,  it  neither  dominates  nor  is  dominated  by  MCF2K,  which  is 
also  weakly  A/T’-hard  to  compute.  As  for  the  SP  relaxation  with  elementary 
routes,  which  is  strongly  A/’-P-hard  to  compute,  it  dominates  all  of  the 
formulations.  This  is  in  accordance  with  Theorem  5.  Finally,  including 
inequalities  makes  a considerable  difference  in  some  cases. 

To  aid  the  reader,  Figure  1 shows  all  known  dominance  re 
all  relaxations  considered.  An  arrow  from  one  relaxation  to  ano  indicates 
that  the  latter  is  stronger  than  the  former. 

6.  Conclusion 

In  this  paper,  we  have  surveyed  the  kno 
ditional  commodity-flow  variables,  and  ini 


armulations  based  on  ad- 
ic<  )?everal  new  multi-commodity 


flow  formulations  that  are  provably  strong  r than  all  of  them,  in  both  the- 
ory and  practice.  Of  particular  interest  is  the  formulation  that  we  have  called 
MCF2b,  which  is  of  polynomial  size  and  yields  lower  bounds  that  are  significantly 
stronger  than  those  obtained  by  all  other  known  formulations  of  polynomial  size. 

A natural  questiqn  fqr  future  research  is  whether  one  could  devise  a compet- 
itive exact  algorithnly|pr  the  CVRP  based  upon  one  of  our  new  formulations, 
or  some  similar  formulation.  To  do  this,  some  kind  of  decomposition  scheme 
might  be  needed.  Another  interesting  research  topic  would  be  to  completely 
characterise  the  projections  into  ai-space  of  the  feasible  regions  of  the  LP  re- 
laxations of  MCFlb  to  MCFld,  MCF2a  and  MCF2b,  and  possibly  devise  efficient 
separation  routines  for  them,  so  that  the  additional  flow  variables  could  be 


Csg£>aratj< 
avoided, 
lar  vehic 


avoided.  Finally,  one  could  attempt  to  adapt  our  formulations  to  other  simi- 
,r  vehicle  routing  problems,  such  as  the  CVRP  with  pickup-and-delivery  [8], 
the  multi-commodity  TSP  described  in  [17],  or  the  VRP  with  load-dependent 
costs  [19]. 
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Figure  1:  Dominance  between  formulations. 
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Appendix 

For  ease  of  reference,  we  include  here  two  of  our  multi-commodity  flow  for- 
mulations: MCFld,  whicu 
that  involves  only  x an 
formulation  of  polynon 


fk(S~m  = fk(S+(k))  = gk(S-(k))  = gk(6+( 0))  =0  (k  e Vc) 


Formulation  MCFld: 


s.t. 


fk(S+(i))  = fk(S~(i))  = gi(S+(k))  = gi(6~(k))  (k,i  G Fc  : i / k) 


kev0\{i,j } 
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